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ABSTRACT 



> I ■ Context. GRB 050721 was detected by Swift and promptly followed-up, in the X-ray by Swift itself and, in the optical band, by the VLT 
03 operated, for the first time, in rapid response mode starting observations about 25 m after the burst. A multiwavelength monitoring campaign 
was performed in order to study its afterglow's behavior. 

Aims. We present the analysis of the early and late afterglow emission in both the X-ray and optical bands, as observed by Swift, a robotic 
telescope, and the VLT. We compare early observations with late afterglow observations obtained with Swift and the VLT in different bands in 
order to constrain the density of the medium in which the fireball is expanding. 

Methods. We have analyzed both the X-ray and the optical light curves and compared the spectral energy distribution of the afterglow at two 
different epochs. 

Results. We observed an intense rebrightening in the optical band at about one day after the burst which was not seen in the X-ray band. This 
is the first observation of a GRB afterglow in which a rebrightening is observed in the optical but not in the X-ray band. The lack of detection 
in X-ray of such a strong rebrightening at lower energies may be described with a variable external density profile. In such a scenario, the 
combined X-ray and optical observations allow us to derive that matter located at ~ 10 17 cm from the burst is about a factor of 10 higher than 
in the inner region. 



Key words, gamma ray: bursts - gamma ray: individual GRB 050721 



1. Introduction 

The Swift Gamma-ray Burst Explorer (Gehrel s"et al. 20041 is 
currently detecting 2-3 gamma-ray bursts (GRBs) per week, 
distributing coordinates with very small uncertainties (few ar- 

Send offprint requests to: L. Angelo Antonelli - e-mail: 
a. antonelli@oa-roma. inaf . it 



cmin down to several arcsec) with delays ranging from few sec- 
onds to tens of seconds after the GRB event. Thanks to its fast- 
pointing capabilities, Swift is able to perform observations of 
the GRB early afterglow phases, both in the X-ray and ultravi- 
olet/optical bands. Moreover, the prompt Swift alerts also allow 
follow-up of GRBs with ground-based facilities. In particular 
the European Southern Observatory (ESO) made its four Very 
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Large Telescope (VLT) units able to quickly react to transient 
sources, allowing them to repoint and start observations within 
just 8 min after the trigger, by developing the Rapid Response 
Mode (RRM) procedure. Before Swift the time needed to deter- 
mine the GRB position was large and most afterglow measure- 
ments could start hours after the burst. Thanks to Swift we can 
now investigate the characteristics of the very early stages of 
the afterglow, when the physical properties of the fireball and 
of the circumburst medium can be derived from the properties 
of the light curve. Moreover, early observations easily span a 
long dynamical range in the afterglow lifetime, so that a rich 
phenomenology can be observed. 

Here we report our analysis of GRB 050721, a long, weak 
burst discovered and located by the Burst Alert Telescope 
(BAT; |Barthelmy et al. 2005) on board Swift and immediately 
(186 s after the burst) followed-up with the Swift narrow field 
instruments: the X-ray Telescope (XRT; Burrows et al. 2 005 ) 
and the UV Optical Telescope (UVOT, |Roming et al. 2005) . 
A previously unknown fading X-ray source was detected 
within the BAT error circle, while no evident optical coun- 
terparts were seen in the optical images obtained by UVOT 
JAntonelli et al. 2005> . The fast accurate localization allowed 
a 0.3 m robotic telescope, located in New Mexico, to de- 
tect the optical afterglow (OA) ~ 369 s after the trigger 
JTorii 2005> . at the level of I c ~ 15.6 mag. The field was 
also observed by the MISTICI collaboration with the ESO-VLT 
UT2 telescope operated for the first time in rapid response 
mode (RRM). The VLT observations started 25 min after the 
GRB and confirmed the presence of the OA within the XRT 
error circle ( Covino et al. 20051 . Such a prompt identification 
of the OA allowed a very dense sampling of its light curve 
at early times, making it one of best examples ever obtained 
(Test a et al. 20 05 1. A spectroscopic observation was also per- 
formed at VLT but the OA spectrum was heavily contaminated 
by the contribution of a bright (R = 16.7 mag) foreground star 
lying very close (1.4") to the OA. No useful information could 
be extracted from it. 

The OA was extensively observed for several days after the 
burst with both Swift and the VLT, obtaining a very good mul- 
tiwavelength coverage. 

In the next sections a detailed study of the XRT and the 
VLT data is presented and discussed. 

2. GRB 050721 

GRB 050721 was detected by BAT on 2005, July 21 at 
04:29:14.3 UT and the BAT on-board calculated location 
was RA = 16 h 53 m 47 s , Dec = -28°23'22" (J2000), 
with an uncertainty of 3' error radius (90% containment, 
lAntonelli et al. 20051 . The refined BAT ground position was 
RA = 16 h 53 m 48 s .5, Dec = -28°23'10" (J2000), with an er- 
ror radius of 3' (90% containment, statistical and systematic). 
The partial coding was 13% ( Feni more et al. 2005 1. 

The masked-weighted light curve shows a FRED-like 
structure with a single large peak starting to rise at To — 5 s 
(Tq being the trigger time), peaking at Tq + 3.7 s and decaying 
back to background levels by Tq + 50 s. The peak is visible in 
the 15-100 keV energy band, but not at higher energies. The 
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Fig. 1. X-ray light curve of GRB 050721 afterglow in the 0.2- 
10 keV energy band. Empty squares and filled circles indi- 
cate WT and PC data, respectively. The curve is background- 
subtracted and the time is referred to the BAT trigger, 2005 Jul 
21 at 04:29:14.28 UT. The solid line shows the best-fit broken 
power law (see Sect. 3.2). The last point is a 3-<x upper limit. 

calculated Tgo (15-350 keV) is (39 + 2) s (estimated error in- 
cluding systematics). 

The power-law photon index of the time-averaged spec- 
trum was 1.81 + 0.08. The fluence in the 15-350 keV band 
was (5.05 + 0.25) x 10~ 6 erg cnT 2 . The 1-s peak photon 
flux measured from To + 3.7 s in the 15-350 keV band was 
(3.4 + 0.8) ph cm" 2 s _1 . All the errors are quoted at the 90% 
confidence level ( Fenim ore et al. 2 005 1. 

3. XRT data analysis 

The XRT observations of GRB 050721 started on July 21, 2005 
at 04:32:20 UT, only 186 s after the trigger, and ended on July 
28 at 08:52:59 UT, thus summing up a total net exposure (in 
photon counting (PC) mode) of ~ 35 000 s spread over a ~ 4 d 
baseline (see Fig. f0- The monitoring was organized in 7 ob- 
servations. The first one was performed as an automated target 
(AT) with XRT in auto state. Therefore, during this observation 
the automated mode switching made XRT take an initial 2.5 s 
image (IM at T Q + 186 s) followed by (T + 194 s) a series of 
windowed timing (WT) frames which were taken until the on- 
board measured count rate was low enough for XRT to switch 
to PC mode (T + 324 s). After this, XRT remained in PC mode. 
The log of Swift observations used for this work is summarized 
in table [T] 

The XRT data were processed with the task xrtpipeline 
(vO.9.9), applying standard calibration and filtering and screen- 
ing criteria, i.e., we cut out temporal intervals during which the 
CCD temperature was higher than -47 °C, and we removed hot 
and flickering pixels. An on-board event threshold of ~0.2 keV 
was also applied to the central pixel, which has been proven 
to reduce most of the background due to either the bright Earth 
limb or the CCD dark current (which depends on the CCD tem- 
perature). For our analysis we further selected XRT grades 0- 
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Sequence 


Obs. mode 


Start time 


End time 


Exposure 


Time since trigger 






yyyy-mm-dd hh:mm:ss 


yyyy-mm-dd hh:mm:ss 


s 


s 


00146970000 


XRT/IM 


2005-07-21 04:32:20 


2005-07-21 04:32:22 


2.5 


186 


00146970000 


XRT/WT 


2005-07-21 04:32:28 


2005-07-21 05:48:30 


205 


194 


00146970000 


XRT/PC 


2005-07-21 04:34:38 


2005-07-21 05:50:49 


1291 


324 


00146970001 


XRT/PC 


2005-07-21 09:00:03 


2005-07-21 19:01:44 


8818 


16249 


00146970003 


XRT/PC 


2005-07-22 00:05:07 


2005-07-22 23:42:59 


19169 


70554 


00146970004 


XRT/PC 


2005-07-23 11:27:27 


2005-07-24 13:13:58 


5937 


197894 


00146970005 


XRT/PC 


2005-08-28 14:50:42 


2005-08-29 08:52:59 


8600 


3320489 



12 and 0-2 for PC and WT data, respectively (according to 
Swift nomenclature). 

3.1. Spatial analysis 

An uncatalogued, fading X-ray source was detected within the 
BAT error circle. We used the xrtcentroid task (vO.2.7) to 
derive the afterglow position, obtaining RA = 16 h 53 m 44 s .62, 
Dec = -28°22'52'.'l (J2000). We estimate the uncertainty to be 
3'.' 3 (90% confidence level). This position takes into account 
the correction for the misalignment between the telescope and 
the satellite optical axis (Mor etti et al. 20051 . 

3.2. Temporal analysis 

During the first orbit of the XRT observation the intensity of 
the afterglow was high enough to cause pile-up in the PC mode 
data. To account for this effect we extracted the source events 
in an annulus with a 30-pixel outer radius (~ 71") and a 4- 
pixel inner radius. For the PC data collected after the first orbit, 
the entire circular region (30 pixel radius) was used, instead. 
The WT data were extracted in a rectangular region 40 x 20 
pixels along the image strip. The selected extraction regions 
correspond to ~ 43% (piled-up PC), ~ 95% (non piled-up 
PC), and ~ 95% (WT) of the XRT PSF. To account for the 
background, data were also extracted in PC mode within an 
annular region (inner and outer radii of 50 and 100 pixels) cen- 
tered on the source, and in WT mode within a rectangular box 
(40 x 20 pixels) far from background sources. Figure^shows 
the background-subtracted light curve extracted in the 0.2-10 
keV energy band, with the BAT trigger as origin of time. The 
last point is a 3-cr upper limit. A fit with a broken power law 
F(t) = Kr ai for t < t h and F{t) = Kt b a ' (t/t b )- a2 for t > t b , 
where % is the time of the break, yields slopes a\ = 2.37 + 0.24 
and a 2 = (1.20 ± 0.04), and a break at t = 399+^ s after the 
BAT trigger Otf ed = 1.39; 52 d.o.f.). The late afterglow (i.e the 
last 3 observations reported in Tabled light curve is fitted by 
a simple power law F(t) = Kt~ a with a = 1 .20*5 ?? ■ 
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Fig. 2. X-ray spectrum of the afterglow of GRB 050721. Top: 
spectra of WT (filled circles) and PC data (empty squares) fitted 
with an absorbed power law model. Bottom: residuals from the 
simultaneous power-law fit to all the data. 
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Fig. 3. Contour plot of the column density vs. the power law 
photon index for the X-ray spectrum of GRB 050721. The con- 
tours refer to 1-, 2- and 3-cr confidence levels. The Galactic 
column density ( 1.60 x 10 21 cirT 2 ; Kalberla et al. 2005 1 is re- 
ported with an uncertainty of 15% (solid vertical lines). 



3.3. Spectral analysis 

Spectra of the source and background were extracted in both 
the WT mode and PC mode in the same regions described 
above at the time corresponding to the first orbit. Ancillary 
response files were generated with the task xrtmkarf within 
FTOOLS v6.0.4, and account for different extraction regions 



and PSF corrections. We used the latest spectral redistribu- 
tion matrices (RMF, v008). The 0.3-10 keV band spectra were 
rebinned with a minimum of 20 counts per energy bin to al- 
low the use of x 1 and fitted by adopting a simple model of an 
absorbed power law model, with the Hydrogen column den- 
sity (at z = 0) kept as a free parameter. In the early part of 
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the observation (first 324 s) the WT mode spectrum yields 
a photon index F = 1.71^^ and a column density Nu = 
2.92+Q 4g x 10 21 cm -2 (slightly higher than the Galactic value: 
1.57 x 10 21 cnr 2 ;|Dickey & Lockman 1990| 1.60 x 10 21 cm 2 : 
IKalberla et al. 20051 . v^ A = 1.00: 55 d.o.f. The observed count 
rate was (10.0 + 0.3) cts s corresponding to an unabsorbed 
flux F x (0.2-10 keV)=7.2 x 10~ I0 ergs s" 1 cm" 2 . In the remain- 
ing part of the first orbit (324-1450 s) the PC mode spec- 
trum yelds a photon index F = 1.86*2}] and a column density 
N H = 3.25+°™ x 10 21 cm- 2 ,xl d = 0.67; 36 d.o.f. The observed 
count rate was (0.72 + 0.02) cts s _1 corresponding to an unab- 
sorbed flux of F x (0.2-10 keV)=1.40 x 10- 1() ergs s _1 cm" 2 . In 
order to enhance statistics, we also performed a simultaneous 
fit of the WT data and the PC data relevant to the first orbit in 
the 0.3-10 keV band. A free constant factor was introduced to 
take into account both the decrease of the mean flux between 
the WT and PC data and the different size of extraction regions. 
We obtain a photon index F = 1.77+gJi and a column density 
Nn = 3.06^ x 10 21 cmr 2 , x 2 KA = 0.86; 93 d.o.f. (Fig.Q. 
Figure shows the contour plot of the column density versus 
the power law photon index, which shows that the column den- 
sity is higher (> 3-cr confidence level) than the Galactic value. 
The PC mode spectrum, obtained by the sum of the last 4 obser- 
vations (see Tab[0, yields a best fit value of F — 2.18^ 22 and a 
column density consistent with the Galactic value (x 2 ed = 0.88; 
18 d.o.f.). The photon index of the late afterglow spectrum, 
compared with the early afterglow spectrum, observed in the 
PC mode, shows a possible softening of the photon index by 
AF = 0.32 + 0.30. 




Fig. 4. VLT image of the field of the afterglow of GRB 05072 1 , 
showing the proximity with the USNO star and the fading be- 
havior of the optical transient. 



4.2. Prompt and late VLT observations 

VLT optical data were obtained in five series during four nights, 
as specified in Table [3] The earliest data came from the acti- 
vation of the rapid response mode, and consist of 33 images 
in the R filter, starting on 2005 Jul 21, 4:54 UT (25 min after 
the GRB) and lasting 53 min. Since the RRM request also in- 
cluded a spectrum, we used the /?-band acquisition image for 
photometry too, hence reaching a total of 34 RRM images in 
the R filter. One hour later, a second series of 4 images in R and 
2 in B were acquired. In the two days after the burst, target- 
of-opportunity (ToO) observations were activated to extend the 
light-curve sampling. In total, the optical sample consists of 39 
measurements in R, 2 in /, and 2 in B. 



4. Optical observations 

4.1. Prompt observation 

The remotely controlled 0.30m telescope, located in the 
New Mexico Skies observatory and operated from the Osaka 
University, observed the field of GRB 050721 starting on 2005 
July 21 at 04:35:23 UT (369 s after the burst), collecting 6 I c 
frames. A visual comparison of the frames with the DSS and 
2MASS images revealed a bright, fading object within the XRT 
error circle. Due to the large pixel scale (1'.'41 per pixel, for 
a field of view of 24' x 24') the optical afterglow was indis- 
tinguishable from the USNO star (U06 16-0429 150) at RA = 
16 h 53 m 44^496, Dec = -28°22'52'.'75 (J2000), (heareafter star 
A, see Fig.|^l. The contribution from star A was estimated from 
a second observation performed on 2005 July 25 with the same 
instrumental configuration and assuming a negligible contribu- 
tion from the afterglow at that time. 

Here we report the first observation only in which the con- 
tribution of star A to the afterglow flux was negligible. The 
observed OA flux was 7 t ,= 15.6 ± 0.2 mag. The field was cali- 
brated using USNO-B1.0 / magnitudes and the total systematic 
error was estimated in absolute photometry as 0. 15 mag which 
is included in the error quoted above. 
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1000 10000 
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Fig. 5. The /?-band light curve of the early afterglow of 
GRB 050721 observed with VLT equipped with FORS2. 
Optical data can be fitted by a simple power law model (F(t) <x 
r") with a — 1.29 + 0.06 (dotted line). This slope is in good 
agreement with the X-ray slope at the same time. 
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Run 


Night 


UT start 


UTend 


Filters 


Images 


Exp. Time [s] 


Seeing 


Observing Mode 


1 


2005 Jul 21 


04:35:23 


05:32:31 


h 


6 


120 


< 1.4" 


Robotic 


2 


2005 Jul 21 


04:53:58 


05:45:21 


R 


34 


30 


0.9" 


VLT+RRM 


3 


2005 Jul 21 


06:40:45 


06:50:35 


B,R 


2,4 


120, 60 


1.2" 


VLT+ToO 


4 


2005 Jul 22 


00:47:11 


04:39:39 


B, R, I 


15, 20, 10 


60, 10 


0.6" 


VLT+ToO 


5 


2005 Jul 23 


02:00:42 


02:07:41 


R 


5 


60 


1.5" 


VLT+ToO 


6 


2005 Jul 24 


00:31:57 


01:51:58 


R,l 


10, 20 


90 


0.7" 


VLT+ToO 


7 


2005 Oct 06 


23:49:37 


00:18:14 


R 


30 


900 


0.65" 


VLT+ToO 



Table 3. Log of R band observations. 



te since burst 


Exp. time 


Magnitude 


Error 


s 


s 


mag 


mag 


1484 


30 


17.93 


0.075 


1559 


30 


18.00 


0.075 


1634 


30 


18.08 


0.075 


1710 


30 


18.17 


0.074 


1794 


30 


18.21 


0.074 


1869 


30 


18.26 


0.075 


1945 


30 


18.33 


0.075 


2021 


30 


18.39 


0.073 


2105 


30 


18.45 


0.074 


2181 


30 


18.50 


0.074 


2256 


30 


18.55 


0.074 


2332 


30 


18.62 


0.074 


2415 


30 


18.65 


0.075 


2491 


30 


18.70 


0.074 


2566 


30 


18.74 


0.074 


2642 


30 


18.81 


0.075 


2726 


30 


18.94 


0.075 


2802 


30 


18.86 


0.074 


2878 


30 


18.92 


0.073 


2954 


30 


18.95 


0.074 


3037 


30 


18.98 


0.075 


3113 


30 


19.01 


0.075 


3189 


30 


19.05 


0.075 


3265 


30 


19.03 


0.075 


3348 


30 


19.08 


0.075 


3424 


30 


19.12 


0.077 


3501 


30 


19.14 


0.076 


3577 


30 


19.20 


0.075 


3661 


30 


19.18 


0.074 


3737 


30 


19.19 


0.074 


3813 


30 


19.29 


0.075 


3889 


30 


19.35 


0.077 


3973 


30 


19.44 


0.077 


4537 


30 


19.48 


0.043 


7891 


240 


20.28 


0.077 


73076 


1000 


21.58 


0.016 


78075 


20 


21.66 


0.062 


85946 


1000 


21.60 


0.024 


248595 


1200 


22.51 


0.110 


6721500 


900 


>25.8 


(5o-u.l.) 



4.2.1. Optical data reduction and analysis 

The optical afterglow was detected at RA=16 h 53 m 44753, Dec 
=-28°22'51'/8 (J2000), very close (174) to a relatively bright 
star (fl=17.4) present in the USNOB1.0 catalog (U0616- 



0429150), see Fig. [4}. In order to avoid any possible contri- 
bution from the nearby star, a PSF profile fitting was adopted 
using DAOPHOT/ALLSTAR dStetson 1987llStetson 1994b within 
IRAF 1 . In order to maximize the efficiency in the detection, all 
the /?-band images from the RRM run were carefully aligned 
and averaged together. Since the FWHM varied slightly during 
the acquisition of the sequence, only the best-seeing frames 
were used to obtain the average master image. The mas- 
ter was then used to create, with daofind, a list of candi- 
date objects, that was used as input list to process the indi- 
vidual images. Thanks to the good mean seeing conditions, 
DAOPHOT/ALLSTAR was able to resolve the two components of 
the pair USNO star/afterglow, providing an excellent PSF fit- 
ting. The output photometry lists were then carefully checked 
for relative zeropoint differences, finding a very good agree- 
ment among the measurements without the need of further re- 
adjustment. All the catalogs from the individual images were 
translated to a homogeneous coordinate system by transform- 
ing pixel coordinates into sky coordinates. The same approach 
was also used later for the images obtained in the other ToO 
series. Magnitudes and, where applicable, colors have been ob- 
tained for all the objects in the master candidate list. FORS 1 
and FORS 2 data have been later matched by transforming all 
the image coordinates to sky coordinates, obtaining a homoge- 
neous catalog. Image calibration was obtained after applying 
aperture corrections to the measured objects; calibration re- 
lations derived from standard star observations were applied. 
Standard stars were observed all nights and the calibrations 
were performed with the IRAF package photcal. Color terms 
have been fitted against the available magnitudes, i.e. B - R for 
the B and R bands, and R—I for the / band, and were found to be 
very small (+0.015 for B-R, < 0.005 for/?-/). Since the stan- 
dards (Landolt fields PG 1323-086 and SA 1 10) were observed 
only once per night, we used the mean extinction coefficients 
for Paranal. 

The VLT prompt observations, obtained in RRM, provided 
a very well sampled light curve. The observed light curve 
shows the typical fading behavior for GRB afterglows well 
described by a simple power law (F(t) = Kr a ) with a slope 
a — 1.29 + 0.06 which is in good agreement with the X-ray 
slope during the same time interval (see Fig.[6ji. A late point, 
obtained at the end of the night (~ 8000 s after the burst and 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 
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~ 6500 s after the first point), is well aligned along the early 
curve. The optical light curve is smoothly decaying and it does 
not show any significant flaring activity. GRB 050721 was also 
observed during the second night starting about 20 hr after the 
burst. We collected our data set in three different images in 
which the afterglow is still clearly detected and we found it 
about 1 .8 mag brighter than the value predicted by extrapolat- 
ing the first night data set. Such a rebrightening was also ob- 
served in the following two observations, obtained respectively 
45 and 69 hr after the burst. Moreover if we consider these two 
points only we obtain a power law decay with a steeper slope 
{a - 1.9 + 0.7) but still consistent with the previous one given 
the large error. A late deep observation (2.5 months later) did 
not reveal either an afterglow or a host galaxy down to a limit- 
ing magnitude of R > 25.8 mag (5cr). We can conclude that the 
observed flattening was not due to the emergence of the host 
galaxy, but that a rebrightening was present in the afterglow. 
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Fig. 6. Multiwavelength light curve of the afterglow of GRB 
050721: from top to bottom / band (black filled triangles), R 
band (red filled squares), B band (blue filled circles), X-flux at 
1 keV (black open circles). Vertical dashed lines show the times 
of the SEDs in Fig. 



5. Results and discussion 

GRB 050721 was promptly observed in both the X-ray and 
optical bands and it was accurately monitored for several 
days after the burst. Such a long term multiwavelength mon- 
itoring showed a different behavior in the X-ray and optical 
light curves. Thanks to Swift, it is now well known that the 
early X-ray light curves of a large fraction of the bursts are 
characterized by a steep decay, followed by a shallower phase 
and finally a somewhat steeper decay (Tagliaferr i et al. 2005| 
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Fig. 7. Spectral energy distribution of the afterglow of GRB 
050721 before the optical bump (f = 8420 s) and at the time of 
the optical bump (f = 76 160 s). Early optical and X-ray fluxes 
are artificially shifted up by a factor of 10 to distance them from 
the later time points. 



INousek et al. 20061 IQ'Brien et al. 2006> . Many events have 
superimposed X-ray flares, probably still related to the 
central engine ( Bu rrows et al. 20051 IFalcone et al. 2 006 
IRomano et al. 20061 |Zhang et al. 2006| IGuetta et al. 2006> . In 
the case of GRB 050721 we did not detect flares in the X-ray 
light curve, while we saw the early steep decay, although 
not one of the steepest so far seen by XRT, followed by a 
shallower phase at about 400 s after the burst. We also did not 
see a second break in the 0.2-10.0 keV light curve, which in 
fact shows a regular fading behaviour up to a time interval of 
about one week, that can be fitted with a simple power law 
(F(i) = Kr a ) yielding a slope of a = 1.20 + 0.04 on the entire 
observation. The optical light curve shows the same regular 
fading behavior over a time interval of about 7000 s with the 
same power law (a - 1 .29 + 0.06) then, at about 24 hours after 
the burst, a rebrightening of about 1.8 mag is observed. 

The consistency between the early optical and X-ray de- 
cay indices suggests that both emissions result from the same 
component during the first hours after the GRB. This fact 
can be futher confirmed by looking at the broadband spec- 
tral energy distribution (SED). Figure shows the SED at 
two times, chosen during the initial decay phase and during 
the optical bump. At the first epoch, the optical data agree 
well with the extrapolation of the X-ray spectrum. The opti- 
cal color {B - R = 1.25 mag) also provides a spectral index 
fiopt — 1.16 ± 0.35, after correcting for the significant Galactic 
extinction {As = 1.2 mag, A« = 0.75 mag; Schlegel et al. 
1998). This value, given its large uncertainty, is consistent with 
that observed in the X-ray range (/i x = 0.86 + 0.17), and with 
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the broad-band optical-to-X-ray slope (J3 x = 0.80). The pic- 
ture is therefore consistent with the optical and X-ray emission 
lying on the same segment of a power-law SED, with slope 
P - 0.80. Thus, no synchrotron breaks were present between 
the two bands. Moreover, the extinction affecting the optical 
data was likely small. 

The spectral and temporal properties of the afterglow 
of GRB 050721 during the first hours are consistent with 
a fireball expanding inside a uniform external medium 
(Sari , Piran, & Narayan, 1998) , providing that the cooling fre- 
quency lies blueward of the X-ray band. For a spectral index 
P ~ 0.80, the predicted decay slope is a — 3/5/2 ~ 1.2, in very 
good agreement with the observed values both in the optical 
and in the X-ray band (their average being a — 1.23 + 0.03). 
The inferred electron distribution index is computed as p = 
2/3+1 ~ 2.6, or, alternatively, p = 4a/3 + 1 = 2.64 ± 0.04. 
Such value is not uncommon among GRB afterglows. 

The situation changes substantially starting ~ 1 d after the 
GRB. At this time, the optical flux is significantly brighter than 
predicted by the early-time decay and the optical spectrum is 
softer (J3 opt = 1.85 + 0.11) than observed at earlier time, and 
steeper than usually observed in GRB afterglows. Both the flux 
and SED suggest that a new mechanism is powering the optical 
emission, leaving unaffected the X-ray region. 

Among the possible explanations for a bump in the opti- 
cal light curve, we consider the emergence of an additional 
energy source like a supernova (SN). It is now well estab- 
lished that long duration GRBs (or a significant fraction of 
them) are produced in SN explosions (Gal ama et al. 1 998 
IStanek et al. 20031 |Hjorth et al. 2003] IMalesani et al. 20041 
Sometime bumps have been detected during the afterglow 
decay days/weeks past the GRB (e.g. IB loom et al. 1 999 
IDella Valle et al. 200 6 ). However, the bump associated with 
GRB 050721 rose very quickly, only ~ 1 d after the gamma 
event, which is much shorter than the observed rising time 
of SNe-Ibc (10-20 days, e.g. |Hamuy, 2003) . One possibility 
is that the SN has occurred before the GRB, as advocated in 
the Supranova model ( Vietri & Stella 1999). An analysis of the 
present data suggest that SNe and GRBs go off simultaneously 
( Delia Valle et al. 2005 ), nevertheless a few days of delay be- 
tween SN and GRB events can be yet accommodated within 
the uncertainties. The major caveat for this hypothesis is rep- 
resented by the magnitude of the bump (R « 21). For a SN 
as bright as SN 1998bw (M R = -19.4; Gala ma et al.l998l . 
the inferred redshift would be z ~ 0.25. This would be diffi- 
cult to reconcile with the nondetection of any host galaxy (our 
limit would correspond to 0.003L* at this redshift). A fainter 
SN would imply a distance even smaller. We thus consider this 
hypothesis unlikely. 

A more exotic possibility calls for the mini-SN scenario 
( |Li & Paczynski 19 98 ), which evolutionary timescale is com- 
parable to that of the rebrightening (1-2 d). However this sce- 
nario would suffer a problem similar to the one faced by a con- 
ventional SN, since the peak magnitude is expected to be com- 
parable to that of a SN. 

A plausible explanation for the rebrightening could 
be given by modifications in the afterglow physics: re- 
brightenings due to energy injection (e.g. refreshed shock, 



IGuetta et al. 20 06 ) or fluctuations in the external medium. In 
the refreshed shock scenario flux fluctuations in both the X- 
ray and optical bands are expected to be proportional to energy 
fluctuations of the fireball. Therefore, a bump in both the opti- 
cal and X-ray light curves is expected. The lack of a rebright- 
ening in the X-ray band rules out this hypothesis. 

In the case of fluctuations in the external medium, we ex- 
pect a bump in the optical light curve (but not in the X-ray one) 
when the fireball encounters regions with enhanced density 
n (e.g. Lazzat fet al. 20021 [Berger et al., 2000). Such a density 
jumps may be caused by winds termination shock, as proposed 
for the case of GRB 030329 ( |Huang et al. 2006) . Nakar, Granot 
& Piran (2003) have shown that the afterglow flux F v at differ- 
ent frequencies may be written in terms of the external density 
in a very simple way. Assuming the typical afterglow parame- 
ters, the synchrotron cooling frequency at 1 d from the burst is 
expected to be v c ~ 10 15 Hz (just between optical and X-ray 
band). Our data on the late X-ray afterglows of GRB 050721 
cannot exclude (within uncertainties) both a steepening in the 
light curve {a — 1 .20^9 ??) and a softening of the photon index 
(AF = 0.32 ±0.30), that are compatible with the hypothesis that 
v c is shifted to energies lower than the X-ray band at this epoch. 
In such a case the X-ray flux (v > v c ) is weakly dependent on 
the external medium density n, while the optical flux (v < v c 
) is F v oc n 3 / 4 (N akar et al. 2 003). Since we see a rebrighten- 
ing only in the optical band (v < v c ) and not in the X-ray band 
(v > v c ) we can conclude that this is an evidence of fluctuations 
in the external medium. In particular we can derive an estimate 
of the variation in the density from the fluctuations in the flux 
m/m oc (F®/Ff ) 4 ' 3 ~ 10. 

However, the optical flux at the second epoch shows a red- 
der spectral distribution than what expected by adopting a sim- 
ple synchrotron model with parameters extrapolated by the X- 
ray flux. A possible explanation for such a red spectrum is 
that the high-density region could be composed of very dense 
clumps, with a covering factor / <k 1 . In such a scenario, the 
rebrightening observed in the optical band may be due to the 
interaction of the fireball with the clumps. Dust in the clumps 
may survive the intense radiation field, so that the spectrum 
produced inside them emerges reddened. Since / <K lthe opti- 
cal flux did not suffer significant extinction, since only a frac- 
tion / was intercepted by the high-density clouds. 

6. Conclusions 

We have presented a comprehensive multiwavelength study of 
the afterglow of GRB 050721 based on RRM and ToO obser- 
vation with the VLT, early robotic observations and Swift XRT 
data. Both the X-ray and optical light-curve are very well sam- 
pled and showed a regular fading behavior (with almost the 
same decay index) within 24 hrs after the burst. At about one 
day from the burst a remarkable (1.8 mag) rebrightening is ob- 
served in the optical afterglow, while no variations are observed 
in the X-ray light curve at the same epoch. A late optical ob- 
servation (2.5 months after the burst) did not show any after- 
glow or host galaxy down to a limiting magnitude of R >25.8 
(5<t) confirming that the observed rebrightening was not due 
to the presence of a bright host galaxy. The broad band SED 
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obtained at two different epochs, during the initial decay and 
during the rebrightening, suggests that during the second epoch 
a new mechanism is powering the optical emission, leaving un- 
affected the X-ray region. Many different scenarios are consid- 
ered in order to explain such a peculiar rebrightening. Since the 
rebrightening was observed only in the optical band and not in 
the X-ray band we propose that this is a plausible evidence of 
fluctuations in the external medium and we derive an estimate 
of the variation in the density from the fluctuations in the flux 
«2/«i ~ 10. In order to take into account an optical spectral 
distribution redder than what expected by adopting a simple 
synchrotron model, we suggest that the observed rebrightening 
may be due to the interaction of the fireball with a high-density 
region composed by very dense clumps, with a covering factor 
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